Quantum effects in a superconducting glass model 



T.K. Kopec* and J.V.Jose 

Department of Physics, Northeastern University, Boston, Massachusetts 02115 

(February 6, 2008) 

We study disordered Josephson junctions arrays with long-range interaction and charging effects. 
The model consists of two orthogonal sets of positionally disordered A'^ parallel filaments (or wires) 
Josephson coupled at each crossing and in the presence of a homogeneous and transverse mag- 
netic field. The large charging energy (resulting from small self-capacitance of the ultrathin wires) 
introduces important quantum fiuctuations of the superconducting phase within each filament. Po- 
sitional disorder and magnetic field frustration induce spin-glass like ground state, characterized by 
. . . not having long-range order of the phases. The stability of this phase is destroyed for sufficiently 

' large charging energy. We have evaluated the temperature vs charging energy phase diagram by 

extending the methods developed in the theory of infinite-range spin glasses, in the limit of large 
magnetic field. The phase diagram in the different temperature regimes is evaluated by using variety 
' ' of methods, to wit: semiclassical WKB and variational methods, Rayleigh-Schrodinger perturbation 

theory and pseudospin effective Hamiltonians. Possible experimental consequences of these results 
are briefiy discussed. 
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C\l ■ I. INTRODUCTION AND MODEL 

; HAMILTONIAN 

, Josephson junction arrays (JJA) have been the sub- 
j — ' ject of considerable recent research activitycl. Such sys- 
I \ terns consist of superconducting grains embedded in a 
0^ ■ nonsuperconducting host and coupled together by the 
Josephson effect. Recently there has been a surge of both 
theoretical and experimental interest in studying disor- 
dered Josephson-coupled systems in an applied magnetic 
fieldHia. Theoretically, it has been shown that in the 
presence of sufficiently strong magnetic fields the system 
^ , may freeze into a state exhibiting spin-sjass-like type or- 
der among the superconducting grainsQ. From the ex- 
\ perimental viewpoint the interest is motivated by the 
O ■ existence of irreversibility lines in the temperature-field 
' diagram, in virtually all high- Tc superconductors thus 
J> , suggesting the existence of a glassy phaseQ. 

For sufficiently small grains, the behavior of each junc- 
^ . , tion in the JJA modified by the quantum effects which 
■ arise from the small capacitances that lead to large charg- 
ing energies. The competition between phase coherent 
ordering and charging effects in periodic JJA has been 
the subject of a number of studiesEl and it is by now 
well established that for sufficiently large charging en- 
ergy quantum phase fluctuations lead to the complete 
suppression of long-range superconducting ordercl. 

In disordered JJA with small capacitances thermal, 
random and quantum fluctuation will determine the 
physics of the system and an interesting question arises 
regarding the competition between them. The problem 
we would like to address in this paper is then: What is the 
effect of having a competition between the thermal, quan- 
tum and random fluctuations on the long range proper- 
ties of a Josephson junction network. This is a highly 
nontrivial question in general but partial answers can be 
obtained in certain limits. Since there arc only very lim- 
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ited studies on this issue, (see, eg Ref. |T^, the purpose 
of this paper is to investigate these quantum-fluctuation 
effects systematically in a mean-field like approximation, 
to be specified below. We shall study a quantum model 
of a disordered JJA system which, in principle, could 
be realized experimentally and simultaneously allowing 
for a detailed theoretical treatment thus constituting an 
attractive setting for the study of the complicated inter- 
play between quenched disorder, interactions and quan- 
tum fluctuations. 

To be specific, we will study a stack of two sets of 
N mutually perpendicular parallel wires (or filaments) 
Josephson coupled at nodes (see, Fig.|^). Since each hori- 
zontal (vertical) filament of the system is directly coupled 
to every other vertical (horizontal) filament, the number 
of nearest neighbors z in this model is z = A^. This is then 
a realization of a JJA with long-range (infinite-range, in 
the thermodynamic limit A — > oo) interactions which dif- 
fer from the conventional 2D Josephson junctions arrays. 
Furthermore, we assume that the distance between neigh- 
boring parallel wires varies randomly around some aver- 
age value I. Finally, the system is placed in a transverse 
magnetic field B and we shall assume that the Joseph- 
son couplings are sufficiently small so that the induced 
magnetic fields are negligible in comparison with B so 
that the phase gradient along any filament results only 
from the presence of the external magnetic field. The dy- 
namical variables of this system are the superconducting 
phases associated with each wire. The properties of the 
(classical^ model defined above have been studied some 
time agon. Very recently the thermodynamical proper- 
ties of this system (both periodic and positionally disor- 
dered version) JiaiSZE been investigated theoretically and 
exp er iment allyM^Ej . 

An important ingredient in our present considerations 
is that we allow for quantum phase fluctuation within 
a wire assuming the node junctions have a sufficiently 
small self-capacitance C. In this case the charging energy 
becomes a dominant quantity considerably affecting the 
properties of the array in the low temperature regime. 
More precisely, we are interested in the behavior where 
the quantum fluctuations compete with the formation of 
the superconducting glassy phase due to randomness and 
magnetic frustration. 

In the low temperature regime we expect to have 
a glassy phase with randomly frozen superconducting 
phases. To examine the extent to which the system 
is frozen, it is convenienjU|to introduce the Edwards- 
Anderson order parameterlla defined by 

QEA = [{S,)'tU, (1) 

where 

S, = [5f,Sf] = [cos(0O,sin(<^,)], (2) 

while (. . .)t and [. . .]av denote thermodynamic and con- 
figurational averaging, respectively. In a disordered state 
the phases will randomly sample their entire phase space 
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and q = 0, while in a frozen system q and [(Si)T]av = 
indicating the absence of long-range order. The spin 
glass-like grjound state may be destabilized by quantum 
fluctuationsCj. As one varies the strength of the charg- 
ing energy (i.e., increasing the strength of the quantum 
fluctuations) there can be a phase transition at zero tem- 
perature between the glassy phase and the paracoherent 
disordered ground states. This quantum phase transi- 
tion in miidom spin systems has received much attention 
recentlyllj. However, its corresponding nature in quan- 
tum disordered Josephson-coupled systems has not been 
previously studied. The goal of this paper is to exam- 
ine this problem in some detail starting from a quantum 
gauge-glass model defined in Eq. (||) below. 

As discussed in the main body of the paper, we need 
to use different complementary calculational techniques 
to attack the problem, for each one of them is by nature 
approximate and their regions of validity are different. It 
is also important that using only one of these techniques 
by itself may lead to spurious results that can only be 
validated by an independent check. A case in point will 
be the reentrant transition found in the variational cal- 
culation, which is valid in the semiclassical region, and 
which is not found in our low temperature expansion. 

We now describe the main body of the paper. In Sec- 
tion II we define the model Hamiltonian, we perform 
the quenched average over the disorder, followed by a 
functional integral formulation of the mean-field theory 
and the derivation of the saddle-point equations in the 
N —>■ CO limit. This constitutes an exact self-consistency 
condition for the glass order parameters. Section III 
presents the phase diagram for the model preceeded by an 
exhaustive investigation of the self-consistency equations 
via a variety of approaches including: semiclassical WKB 
and variational methods, Rayleigh-Schrodinger perturba- 
tion theory and pseudospin effective Hamiltonians in a 
truncated charge Hilbert space. Finally, in Section IV, 
we discuss our results and present our conclusions. 

II. THE MODEL 

The quantum Hamiltonian of the disordered system of 
Josephson-coupled wires is given by 

Ej ^ ^ 

^■'^ 7^ XI 51 ~ cos(0ft, - (t)yj - f,j)] . (3) 

Here, hi = {2e/i)d/d(t>i is the charge operator while 
{4>vj) represent the superconducting phase of the i— th 
horizontal (i— th vertical) wire, respectively; Ej is the 
Josephson couphng and K = Ae^ /C the charging, respec- 
tively. In order to have a well defined thermodynamic 
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limit we have to scale the Josephson coupling energy by 
a factor VTV- Furthermore, fij = (27r/$o) //.^ A • dl is 
the line integral of the vector potential A and $o is the 
elementary flux quantum. 

It is clear that Ej is a positive quantity (and it may 
depend only on the distance between wires which we are 
neglecting) and it does not introduce any frustration in 
the system. Magnetic field and random location of the 
wires is what generates variations of the phase parame- 
ters fij , thus allowing for the random frustration present 
in the systemLD. The relevant quantity is thus the effec- 
tive random coupling matrix in the Josephson part of the 
Hamiltonian iS) 



Ej 











whose density of the eigenvalues is given by 



p{E) ^ -(l/7r)Im[£;5,, 



(4) 



(5) 



where the bar denotes averaging over the positional dis- 
order. The behavior of varies with increasing the 
strength of the magnetic field. In the case of large 
magnetic field, so that the flux per average plaquette 
(f> — Bp is much larger then the elementary flux quan- 
tum ^/^o S> l,the frustration parameters fij acquire 
random values and fill the interval (0, 2tt] uniformly. In 
this limit the effects of disorder become especially ap- 
parent and the behavior moves toward the asymptotic 
regime which is essentially field independents. In this 
case correlations between the matrix elements van- 
ish and the density of states approaches the Wigner 
semicircular lawn 



p{E) pw{E) = (TrEjy'Jl - iE/2Ej) 



(6) 



implying that the random matrix (|^) belongs to the 
Gaussian unitary ensemble with only the second moment 
non- vanishing, 



(JiiJkl)av — -^SikSjl. 



El 
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(7) 



Therefore, in the high-field regime, we can implement 
a mean-field theory of our quantum gauge-glass prob- 
lem in a way closely resembling the infinite-range in- 
teraction Sherrington-Kirkpatrick magnetic spin-glass 
modellS. However, the resulting formulation is not just 
a quantum extension of the planar spin glass model. Un- 
like the random bond XY spin glass, the improper global 
rotation (j)i —f —(pi (i.e. time reversal) is not a symmetry 
of our quantum gauge-glass model. To make this obser- 
vation apparent we introduce the two component "spin" 
vector (H) so that the Josephson part of the Hamiltonian 
(0) reads 



H, = - 



Ej 



N N 



2VN 



(8) 
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where z is a unit vector perpendicular to the plane con- 
taining Si. We note that, in addition to the conventional 
XY coupling Shi-Svj there is a cross term Shi x Syj which 
is the analogue of the spin-orbit (i.e. Dzialoshinsky- 
Moriya (DM)) interaction in magnetic systems - essen- 
tially violating the time- reversal symmetryllj. Thus, the 
present quantum gauge-glass problem formally resembles 
more closely the quantum spiifcglass formulation in the 
presence of the DM anisotropyEj. 



III. DISORDER AVERAGE AND MEAN-FIELD 
FORMULATION 



In a random system we need to calculate the average 
of the free energy density = —{1/(3) \nZ/2N over the 
disorder (^. This is done by using the replica method 
permitting to average the replicated partition function 
Z" instead of In Z, 



f3T = ~ lim lim -^((Z'-' 



1) 



(9) 



It is convenient to express the replicated partition func- 
tion Z" = Tr exp(— /377") in the interaction repre- 
sentation as 



Z" = Tre~'^S»-f^''T^exp 



f drHfir) 



with the interaction picture Hamiltonian 

Hf{T) = e^"«Hfe-^"^ , 
and the free part 
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(10) 



(11) 



(12) 



where H°' = Hq + Hf is the total Hamiltonian. For the 
interaction Hamiltonian one has explicitly 



Hf{T)^ ^ j E E ^^^'^ [S^.(^) • S.",(r) - *z . S^r) x S^^ir)] + h.c | , (13) 



N N 



so that the statistical average can be taken in the en- 
semble given by H^. Here Tr is the Matsubara "imag- 
inary time" ordering operator allowing us to treat the 
time dependent operators Sf (r) = e~'^^°Sfe'^^° as c- 
numbers within the time-ordered exponential ( p^ ) . Con- 
sequently, the Gaussian average (0) readily gives 



(Z")^^ = Tre-So'3«oT^exp 



dr dT'n{T,T') 

Jo Jo 



(14) 



where 
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i—1 j — 1 0.(3 



-[S^.(r)xSe,(r)].[Sf,(r')xS^^.(r')]}. 



To make further progress we utilize the vector identity 
(A X B) • (C X D) = (A • C)(B • D) - (B • C)(A • D) 

to reduce mixed vector products along with integrations 
over auxiliary variables according to the formulas 



exp(±a6): 
exp(aV2)= 



dxdy 
27ri 
dx 



e^p[±{xy -ax~ by)], 
eyip{~x'^/2 — ax), (16) 



allowing to decouple various quartic interactions in 
Eq.(|l5|) and reduce Eq.(p^) to the effective single- 
filament problem. In terms of the functional integrals 



(^)av = J l[l[DQf^DR'^^DP-^iDP-^iexp{-2NL[I>,,I>2,Q,'R]), 



involving the non-local (in time) tensor fields 
X^^^(r, t') (X = Pi, P2, Q) and i?"^(r, r') where /z, 1/ = 
x,y. The effective local Lagrangian reads 



L[Pi, P2, Q, R] = TrPiPa + TrQ^ + TrR^ - ln$[Pi, P2, Q,R], 



with 



TrXY= f'dr T dr' ^ ^ (t, r')rj;"(r', r), 

-^0 a/3 

TrR2= / dT / dT'y^R''^{T,T')R'^''{T',T). (19) 
Jo Jo „« 



Here, 



$[Pi,P2,Q,R] = Tre-'^^T.exp 

is the effective time-dependent single-filament Hamil- 
tonian describing interactions between replicas 



dT / dT'H,s{T,T') 



HMr,r')^~EjY,J2 



aP fJ.i' 



Q^.^(r,r')-iPif.(r,r')+ 



R'^^{r,T')S,^~lp^,l{T,T') 



St{r)SP{r')- 



(21) 



In the thermodynamic limit, N ^ 00, the steepest de- 
scent method can be used which amounts to finding the 
stationary points X^^^^ and i?g ^ determined by the ex- 
tremal conditions 



5L[Pi,P2,Q,R]/^X;^,^ = 0, 
JL[Pi,P2,Q,R]/^i?"'' = 0. 

Thus, one obtains 



(22) 



QtUr-T')=^GfAr-r'l 

Po,t.ir r')= P,%^{r - r') = ^Cf^r - r') (23) 

where the correlation function describing the dynamic 
self-interaction is 



Tre-/'«or,5^(r)55(r')exp 


[-i 


/3 


f dT2Hcs{Tl,T2) 



Tre~^^T^exp 




('^ dTl 



I dT2Hcs{Tl,T2) 
lO 





(24) 



Equation ( |24|) represents an exact self-consistency con- 
dition for the replica dependent matrix Green function 
G"^ [t — t') of the quantum disordered Josephson model 
(^. For classical spin glasses, this is a matrix G"'', where 
the off-diagonal components of q^p can be related to the 
spin-glass-like Edwards Anderson order parameter 

qEA= maxQ^^(9"^), 

g"/3=G«^(r-r')<5^., (25) 

which is purely static (i.e., "imaginary-time" inde- 
pendent) and vanishes on the superconducting-glass- 
paracoherent phase boundary. However, for the quantum 
problem the time dependent fluctuations of the replica 
diagonal components GJ|J"(t — r') must be considered in 
the " imaginary" Matsubara time r. More precisely, the 
replica diagonal part G"" is not an order parameter be- 
cause its expectation value is non-zero on both sides of 
the transition; nonetheless it has to be determined self- 
consistently along with the glass order parameter g"'^. 



IV. CALCULATION OF THE PHASE DIAGRAM 



A general solution of the self-consistency equation 
( p^ ) poses a rather difficult problem since the quantum- 
mechanical nature of the problem requires that the 
time-dependence of replica-diagonal dynamic parame- 
ters G|^"(r — t') have to be determined self-consistently. 
Therefore, we employed here the static ansatz (cf. Ref. 
^1|) which retains only the uji^q Fourier component 

tS^v = {1/(3) Jq G'/^di''') of the dynamic self-interaction. 
Furthermore, since we are here interested mainly in the 
critical line separating the glass and paracoherent phases 
(where qea = 0) we employ the replica-symmetry as- 
sumption {q"P = q for a ^ /3). In this way we avoid 
the subtle intricacies of the replica symmetry breaking 
(which, however, will be important inside the glass phase 
region). Therefore, 

G';:S{T) = [rS^p + q{l-S^f,)]S^, (26) 
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and the effective Hamiltonian becomes 



i?eff(r, r') = E E [2^^"/^ + 29(1 - 5o.p)] S';:{t)SI{t'). (27) 



2 

a/3 M 

Consequently, the averaged free energy density (^, in 
the rephca n — > hmit, becomes 

Jo 

/>oo 

Jo 

Z{cr, p)= Tr^exp [-/^^^(fJ, p)] , 
' " 2 902 + 

-£;,/V2 (ctV^ + pV^^) cos((/.). (28) 



^^(cr, p) = 



Here we have employed integrations over the auxiliary 
variables cr, p and Tr^ . . . — • ■ ■ Y^ii't')), with 

|\E'£((/))) being the eigenstates of the operator H^{a,p) 
and 



roo 

q 



poo 

r= / 
Jo 



adae-^'/^icos^ {(/))) a (29) 



where 



1 

xTr^...exp[-/3H^((T,p)]. (30) 

The freezing temperature, i.e. the onset of the glassy 
phase, is marked by a non-zero value of the spin glass 
order parameter q. We can now establish the equation 
for the critical line Tc{K) by expanding the free energy 
in powers of q and equating the coefficient of q^ to zero. 
We find that 

P,Ejr{P,,K)^l, (31) 
and the self-consistency equation for r is (cf. Eq.(^8|)) 



pdpe-P /^Tr^ cos2(0) exp [-PU^^P, 0)] 

r='-^^^ (32) 

/ pdpe-p"'^TT^ exp [-m<p{p, 0)] 
Jo 



with the effective single site quantum rotor Hamiltonian 

y 



^0(P,O) = - — — -Aicos(0), (33) 



where p = Ejp\j2r. Finally, using partial integrations 
it is convenient to represent the self-consistency equation 
for r as 
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poo 

/ pdpe-p'/^Zip, 0) [2{l3Ejfr^ + 3 - p'] = Z(0, 0), 

(34) 

where Z(0,0) = 6*3(0, /?/s:/2), and 63(0, PK/2) is the 
theta function 

00 

6'3(z,<) = l + 2^cos(2fcz)e-'='*. (35) 

k=l 

To proceed, we have to exphcitly evaluate Z[p,0) which 
amounts to finding the eigenenergy of the effective quan- 
tum rotor Hamiltonian ( |33[) as a function of (arbitrary 
positive) p. This is a standard eigenvalue problem in- 
volving a Mathieu type differential equation. Unfortu- 
nately, it is difficult to obtain a useful analytical solution 
for general values of K and p in this way. Therefore, we 
treat the problem in various charging energy-temperature 
regimes using a combination of calculational approaches 
to construct the entire phase diagram. 

A. Semiclassical WKB limit 

For ksT ^ Ej, K/Ej <C 1 the charging energy is ex- 
pected to play a subdominant role. Our approach here is 
to consider the "potential energy" p cos{(j)) in the quan- 
tum rotor Hamiltonian ( ^3| ) as the leading contribution 
and treat charging-energy perturbatively. At lowest order 
in K one obtainf£3 

{p) « -EjpV^^ [1 - {(3K/2A)] cos{(f) . (36) 

Furthermore, by using the identity 

InW - /'"^e^™^(*)cos(n0), (37) 
Jo 2tt 

where /„(A) is the modified Bessel function of nth-order. 
Using Eq. ( |34| ) , the critical boundary for small K is read- 
ily obtained from 

/•oo 

/ pdpe-P'-/'{5-p')l, 
Jo 

(38) 

For K = at the classical critical point we have 
ksTc/Ej = V2/2 « 0.71 (cf. Ref. H). 

B. Variational method 

The quantum mechanical partition function Z{p, 0) in 
Eq. ( ^^ can be approximated by an effective classical 
function via the path integral formalism using a varia- 
tional method. Following Refs. [2^:^ (see Appendix for 
details pertaining this problemjwe obtain the critical 



,V2/3Ej ( 1 



(3K 
~2A 



boundary between the glassy and paracoherent phases 
from 



smh[(3Kn4p)/2] '° 



where ^c, ^^'^ determined (for a given value 

of p) from the set of self-consistency equations: 



Up) 

„2 



2 " " 

1 



coth 



^l=P 



■expl-- 



V 2 



(40) 



jexp 



(39) 



Numerical evaluation of Eq.(p9|) reveals a reentrance 
in the low temperature region from the glassy phase 
back to the paracoherent state (see Fig. 2) and the zero- 
temperature critical value of reduced charging energy 
ttc = K/Ej « 1.1 In the opposite semiclassical limit, 
A' — > Eq.(^ reduces to Eq.(H) (see Appendix). 
The important question arises as to whether the pre- 
dicted reentrant feature is a genuine property of the 
model or it is an artifact of the approximatipiif—We note 
that the variational method implementationBjO consid- 
ers that the paths must satisfy the boundary condition 
^(0) = 4>[(3), appropriate for a quantum particle in a 
periodic potential rather than a quantum rotor system 
(Eq.(^3|)) with 27r periodic wave functions. Therefore, 
the proper boundary conditions for the latter must be 
0(0) = (/)(/?) -I- 2nm (where m = 0, ±1,±2,... are the 
winding numbers) which appear not to be accounted for 
consistently in this approach. This might have conse- 
quences for the low temperature behavior of the system 
and we are therefore motivated to look for an alternative 
method to study the quantum T ^ limit. 

As we shall see next, the reentrant transition found 
in the variational approximation is not found in the low 
temperature expansion. Thus it is most likely an arti- 
fact of the variational approach, which is strictly a high 
temperature approximation. 



C. Perturbation expansion about the quantum limit 



Assuming thati^T ^ Ej the potential energy of the 
quantum rotor ( |3^ ) pcos{(j)) may be treated perturba- 
tively using the standard Rayleigh-Schrodinger approach. 
The unperturbed part is then —{K/2)d^/d4)^ with eigen- 
functions 

|C('/'))-^exp(zm0). (41) 
To lowest nontrivial order we get 
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K 

y 



^4 



rp 



1 - AmP- 



(42) 



SO that 



m— +OC 



z(p)= 



-/3£;™(p) 



(43) 



m— — ctt 



and the critical boundary condition is implemented after 
performing the integration obtaining 



m— +00 



l&m^ + 8m2 [i{Ej/K) - 1] - 8{Ej/K)^ - Q{Ej/K) + 1 
[4m2 + 4:{Ej/K) - if 



(44) 



At T = the TO = term is the only one that con- 
tributes and thus 



1 



l-AiEj/K,{T^O)) 



l-A{Ej/KrXT ^Q)) 



= 1, 
(45) 



that gives a^T = 0) = (6 + 2VT7)/2 « 7.1231 which 
differs quite significantly from the critical vaLne-of a ob- 
tained previously by the variational methodllil'tj. We at- 
tribute this difference again to the fact-thai in the present 
formulation of the variational methodcja the boundary 
conditions for the superconducting phase variables are 
not compatible with the discrete nature of the charge 
transfer process which becomes especially important as 
T — !■ 0. As a result, a discrepancy in the low temperature 
limit is expected, whereas in the classical limit [K 0) 
Tc is reproduced correctly. 



D. Truncated charge space projection method 

The expansion about the quantum limit from the pre- 
vious subsection is not sufficient to treat the critical 
boundary away from the quantum critical point adT = 
0), where the condition a S> 1 is no longer valid. A non- 
perturbative approach is sought and of particular interest 
are the truncated charge state models (TSCM) spanned 
by the charge states of the operator —{K/2)d^/d(j)^ in 
a restricted, finite-dimensional, Hilbert space. This can 
be interpreted as an approximation of a suitable quasi- 
spin model. Using the eigenstates of the charge operator 
( [4l| ) one can readily calculate the matrix elements of the 
operators involved 



Nkr,i= ($fe((^)|n((/))|$,„((/))) = J —exp{-ik(j)) f — ^ ) cxp(im0) 2emSkm, 



[S.]km= {M^)\cos{miW) ^ I ^exp[-z(fc-m)0]cos(0) 



2tt 







2('5fe-m-l,0 + Sk-m+lfi). (46) 
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In particular, for k.m = 0, ±1, one has 



Nkm— [Sz]kmt 

where Sa{a = z, x) are given by 



(47) 




5. = |00 0|, 5, = -i=|10l|. (48) 



1 

1 1 
1 



Consequently, the lowest-order quasi-spin model belongs 
to 5 = 1 and its Hilbert space is spanned by the charge 
states |0), I ± 1). It is now useful to recast the quantum 
rotor Hamiltonian (E3) in the 5 = 1 pseudospin language 



H - {Kl2)Sl - (/i/V2)5,. (49) 

In the context of magnetic spin glasses the first charg- 
ing energy term in Eq.(^9[) refers to a single-ion crystel 
anisotropy, which opposes ordering in the x — y planecJ. 
Thus the system will exhibit a phase transition driven 
by quantum fluctuations of the transverse (pseudo)spin 
component. Correspondingly, for the statistical sum 
Z(p, 0) we have 



+2exp(-MVoshK^^-^)^V^-'^^^^^M, (50) 



while 

poo 

pdpe-p"/^Z{p, 0) (5 - p^) - Z(0, 0), (51) 



and the critical line TJa) can be readily calculated nu- 
merically using Eq.(^ll) (see Fig.||). We found ac{T — 
0) = 7 in good agreement with the perturbative ap- 
proach. The classical value Tc{a = 0) is underestimated 
as a result of restricting the original charge state Hilbert 
space. We have also examined models spanned by five 
(|0), I ± 1), I ± 2)) and seven (|0), | ± 1), | ± 2), | ± 3)) 
charge states which better approximate the behavior at 
high temperatures. For both cases it is possible to derive 
critical line equations analogous to Eq.(|5^) analytically. 
However, the corresponding formulae are too lengthy 
to reproduce here. We found that close to Tc{a — 0) 
TCSM exhibit a small reentrance for certain interval of 
a's. However, as the number of charge states mchargo 
increases this interval became narrower and presumably 
disappears for mchargo — > oo indicating that this might 
be a spurious feature due to the restriction imposed on 
the original infinite-dimensional Hilbert space of charge 
states. 
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V. CONCLUSIONS 



In this paper we have studied the competition be- 
tween quantum and thermal fluctuations in a supercon- 
ducting glass state using positionally disordered model of 
ulthrathin Josephson coupled wires in a transverse mag- 
netic field. We focused our attention in the regime where 
the magnetic field produces large frustration in a sys- 
tem of randomly spaced Josephson microj unctions with 
large number of nearest-neighbors z = N . We inves- 
tigated this model by spin-glass inspired techniques in 
the N ^ oo limit, showing that that the model is a 
superconducting analogue of a quantum spin-glass with 
Dzialoshinsky-Moriya time-reversal braking interaction. 

We studied the phase boundary separating the para- 
coherent from glassy phases as a function of the charging 
energy associated with the small capacitance of an indi- 
vidual wire. We found that, for sufficiently large charg- 
ing energy, the glassy ground state is destabilized due to 
the strong quantum fluctuations. This is reminiscent of 
the scenario found in ordered Josephson junction arrays, 
where charging energy effects can lead to the destruc- 
tion of long-range phase coherence by zero-point quan- 
tum fluctuations. However, this analogy is not complete 
as in the glassy phase there is no long-range order and, 
therefore, the destructive role of quantum fluctuations is 
less transparent. We note that it is not easy to improve 
our static ansatz analytically. We expect that a direct 
Quantum Monte Carlo calculation, where the imaginary 
time direction is discretized, will yield information inside 
the phase boundary. The shortcoming of this approach 
is that one can not get too close to the T = region. We 
leave this problem for the future. 

The results presented here can be tested experimen- 
tally in artificially fabricated disordered arrays of ultra- 
thin filaments provided that the charging energy of a 
wire is large enough to produce substantial supercon- 
ducting quantum phase fluctuations. For example, Tc{a) 
should be observable by resistivity measurements. In ob- 
taining our results for the phase diagram we considered 
the solution which does not break the replica symmetry. 
Although a broken replica symmetry solution is not re- 
quired to trace the critical phase boundary (where the 
order parameter vanishes) it is important when distin- 
guishing equilibrium from nonequilibrium properties of 
the glassy phase. In the low-temperature phase, as usual, 
the strongest signatures of glassiness are in the dynami- 
cal properties. Similar to the magnetic spin glasses, non- 
ergodic behavior is likely to manifest itself in differences 
between field-cooled and zero-field prepared samples via 
typical effects like hysteresis, remanence, aging effects 
etc. In the classical limit of our model (vanishing of 
charging energy) all these history jd^etident signatures of 



ing non-equilibrium metastable states might be probed 
eg. via ac conductivity measurements determining the 
barriers separating metastable states and the associated 



a glassy state have been observedtHI 




I. The correspond- 



13 



distribution of relaxation times as a function of charging 
energy. 

Although the system of superconducting ultrathin 
wires constitutes an interesting experimental setting to 
test the predictions of the infinite-range interaction the- 
ory, the investigation of the behavior of disordered quan- 
tum truly 2D array would require knowledge relevant to 
short-range spin glasses and presents a difficult subject 
for further study. 
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APPENDIX A: VARIATIONAL METHOD 



In this appendix we give the specific details of 
the derivatiop of Eq. (39) based on the variational 
approachoE3. The path integral formulation of the par- 
tition function involves an infinite product of periodic 
paths (/)(t) =00-1- J2'eLi{4'iS^'^''^ + c.c) in the form 



Z(p,0) = 



n 



ire J J.im 



7:/iPKLuj)_ 



X exp • 



dr cos 



The essence of the variational method 



(00) 



is to ap- 
proximate (Al) by an effective cZassica/ partition function 

Z(p,0)= /^^e-'3^--(^«). (A2) 

Here, the associated effective potential V^s is given by 



= i In 



sinh(/3ifrj/2) 



PKn/2 



- ^^n^a\ 



(A3) 



The unknown functions a^((/)o), 1^2(00) and il(0o) are 
determined by using extremal principle from the self- 
consistency conditions 



1 



Va2 



l3Kn , /(3Kn\ 



1 



%/2^ 



: exp 



2a2(0o) 



cos{(t>), 



;^exp ( - — ) cos(0o) 



(A4) 



(Al) 
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For low temperatures w e se ek for the uniform solution 



(i.e., with (jjo = 0) of Eq.(|A4D (cf. Eq.(y)). In the regime 



of high temperature and small quantum fluctuations by 
using the expansion 



-coth(-) -l] = —~ — +0(x^), (A5) 
2 V2/ J 12 720 ^ ^ ' 



1 

X 

we have for Vcs{4>o) 



n{p,0)^ -EjpV2rcxp ( cos((/.o), 

EjpV2r[l~^)cos{M, (A6) 



24 , 

i.e., T4ff(0o) reduces to semiclassical WKB result ( 
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FIG. 1. Disordered Josephson-coupled array of 2A'' super- 
conducting wires (straight horizontal and vertical lines, square 
box denotes the junction at a node). Each wire has small 
self-capacitance C and is characterized by the superconduct- 
ing phase (j>i- The magnetic field B is applied perpendicular 
to the array. 
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FIG. 2. Temperature-charging energy, T — a, phase dia- 
gram of the disordered Josephson-coupled array of supercon- 
ducting wires with self-charging energies. Hero a = with 
K the charging energy and Ej the Josepiison coupling energy. 
The results were obtained from variational method (VM) and 
truncated charge states models (TCSM) for different num- 
bers of charge states as indicated. We note, as discussed in 
the text, that the reentrant behavior obtained from the varia- 
tional calculation does not emerge from the low temperature 
analysis and thus it must be an artifact of the variational 
approximation. 
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